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We review the theory of meson decays to two lepton pairs, including the cases of identical as well 
as non-identical leptons, as well as CP-conserving and CP-violating couplings. A complete lowest- 
order calculation of QED radiative corrections to these decays is discussed, and comparisons of 
predicted rates and kinematic distributions between tree-level and one-loop-corrected calculations 
are presented for both vr" and K*^ decays. 
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I. BACKGROUND 

Meson decays to two photons should exhibit interest- 
ing correlations between the photon polarizations Al- 
though existing particle detectors cannot measure pho- 
ton polarizations directly, it has long been known that 
the polarization correlations can be measured indirectly 
by studying angular correlations in the related double 
Dalitz decays in which both photons undergo internal 
conversion to a lepton pair. More recently, it has been 
pointed out§ that a detailed study of these correlations 
can be used to determine the relative amount of two pos- 
sible meson-77 couplings (one CP-conserving and one 
CP-violating for mesons that are CP eigenstates) that 
can contribute to this process. 

Dalitz and double Dalitz decays are also of interest be- 
cause they can be exploited to perform a measurement of 
the electromagnetic form factor of the decaying meson — 
that is, how the meson couples to one real and one vir- 
tual photon (Dalitz decay) or two virtual photons (double 
Dalitz) depends on the values of the photon(s). An 
accurate knowledge of this form factor is essential, for ex- 
ample, to calculate the so-called long-distance contribu- 
tion to the rare decay Kl M"*"/^"- The short-distance 
contribution to this process, mediated by loops involving 
heavy quarks and massive vector bosons, is sensitive to 
the CKM matrix element Vtd, but this contribution can- 
not be extracted from the accurate experimental mea- 
surement of the partial width unless the long-distance 
amplitude is precisely known. 

The tree-level rates for several double Dalitz processes 
have been published in various forms l^], |3[ ^, ||. The 
first experimental observation of a double Dalitz decay 
was published in 19620 |^. A total of 206 examples of 
the decay tt^ — > e'^e^e^e^ were observed by Samios in 
a sample of some 800,000 bubble chamber photographs. 
Based on the observed angular correlations, Samios was 



able to exclude the possibility that the 7r° was a scalar 
particle (with a CP-conserving decay) at the 3.3a confi- 
dence level. His measurement of the branching ratio for 
that process remains the only one published to date. 

Experimental observations of the much rarer kaon dou- 
ble Dalitz decays began to appear in the 1990's. Sev- 
eral measurements have been made of the decay — *■ 
e+e-e+e-|, |, the most recent of 

which are based on several hundred observed events. The 
still rarer double Dalitz mode — > e~^e~ fj,'^ ^~ is par- 
ticularly interesting because it is free of complications 
arising when there are two identical lepton pairs in the 
final state, and because it probes only the kinematic re- 
gion in which one of the virtual photons has > 4to^. 
The first example of this decay was reported in 1996p5[. 
In 2001 the KTeV experiment published a branching ra- 
tio based on a sample of 43 events most recently, 
KTeV has reported results from a combined sample of 
132 events, including the earlier 43|]l7|. 

Further experimental results on both ir^ — > e^e^e^e^ 
and the two kaon decays are expected in the near future 
from the NA48 and KTeV experiments. As the statistics 
available to the experimenters increase, it will be nec- 
essary to have a more accurate theoretical description of 
these decays, incorporating the significant efi^ects of QED 
radiative corrections. These corrections, discussed in this 
paper, have a significant impact on the extraction of both 
form factors and angular correlations from high-statistics 
data. 



II. TREE-LEVEL AMPLITUDES 

The most general CPT invariant interaction governing 
the transition of a spin-zero meson to two photons is 
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(1) 

where J^p and J^s are dimensionless form factors for a 
pseudoscalar and scalar coupling which may be momen- 
tum dependent, F^'^ is the electromagnetic field tensor, 
and 6 is the field of a meson of mass M. The factor 
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of —i defines a phase convention in which the form fac- 
tors will be real if CPT invariance holds and there are 
no absorptive decay amplitudes. It will be convenient to 
decompose the couplings into real and imaginary parts, 
and define them (up to an arbitrary overall phase ^E") in 
terms of a mixing angle C, and a phase difference S, 

Tp = gpfp{kl ki)e^^ = ~gfp{kl kl) cosCe'*, (2a) 
^5 = 55/5(fc?,A:|)e^* =5/5(fc?,fc2')sinCe''e^*, (2b) 

where = gp + g| and (5 = <f> — '5. fp and fs are 
functions of the k"^ values of the two virtual photons and 
are normalized such that /(0,0) = 1 for both couplings. 
In principle it is also possible for 5 to depend on k\ and 
^2- 



A. Two— Photon Decay 

The two-photon partial width provides information 
about the coupling constant g along with some con- 
straints on the mixing angle C and phase difference 5. 
The matrix element for the decay to two real photons 
with helicities Ai and A2 is 



M 



A1A2 



M 



(3) 



where k^ and are the momentum and polarization of a 
photon. The calculation of the two-photon couplings for 
photons of arbitrary mass is carried out in Appendix 
For real photons, the kinematic factors in the couplings 
reduce to A = z = 1 and w = 0, and the momentum de- 
pendent functions fp and fs reduce to unity. Therefore, 
one has the following two contributions 

M++ = -Mg (sin Ce*'' - i cos C) , (4a) 
M— = -Mg{smCe'^ + icosC) . (4b) 

The partial width is obtained by integrating 1/(2M) 
times the squared matrix element over the available phase 
space. The matrix element is a constant, so the integra- 
tion results in a factor of l/(167r). The partial widths for 
the two allowed helicity states are then 



;;2 



Alg 
327r 
Mg^ 
32n 



(1 — 2sinCcosCsin^) , (5a) 
(l + 2sinCcosCsin(5). (5b) 



The decay rates to the two final states will be identical if 
either the phase difference between the two couplings is 
zero or one of the form factors is zero. In any case, the 
total rate is equal to T^j = Mg^ / (IGtt). An experimental 
measure of the two-photon width then gives a value of 



B. Application to Neutral Pion and Kaon Decays 

One can compute g for the 7r° 77, Kl — > 77, 
and Ks 77 decays using the current experimental 
values of the two-photon branching ratios along with 
the meson lifetimes and masses, with the results 



go = (1-70 ±0.06) • 10"^ (7a) 
5L = (8.75 ±0.12) • 10^1°, (7b) 
5s = (1-36 ±0.11) • 10"^ (7c) 

It will be useful to define CP-even and CP-odd 77 states: 



l(77)±) = (l++>±|--»/V2. (8) 
The matrix elements for tt" — ^ 77 are then given by 

((77)+ 1 T Ivr") = -ZN/2Af,gosinCoe**°, (9a) 
((77)- 1 T IttO) = -\/2Af,5oCOsCoe'*". (9b) 

Similar expressions can be given for the Ki and K2 de- 
cay matrix elements with the subscripts changed to 1 
or 2 respectively and M.^ replaced by Mk- These ex- 
pressions may be compared to those in Ref. ]l9| in which 
CP-violation observables for the if ^ 77 decays were 
calculated. In that article Sehgal defines 

((77)+l T \K,) 
((77)-|r|ifi) 
((77)+ 1 T \K2) 
((77)_|r|/f2) 





(10a) 


ICqC ^ , 


(10b) 




(10c) 




(lOd) 



In our notation we have 



Ce = 
Co = 

de = 
do = 



V2M_R-.gi sin Ci, 
~V2MKgi cosCi, 

V2MKg2 sinC2, 
-V2MKg2 C0SC2. 



(11a) 
(lib) 
(He) 
(lid) 

Oe~Po 



(6) 



The observable phase differences are given by 5i = 
and S2 = He — fJ'o- 

The phases 6i are observable only if CP is violated in 
77 decays. As Sehgal noted in Ref. |l^], the phase So in 
7r° decays should be very close to zero if CPT is conserved 
because the relevant absorptive amplitudes are of order 
a^. In the case of kaon decays, on-shell intermediate 
states such as tttt couple strongly so that the phases Si 
and S2 may be large even if CPT is conserved. 



C. Four— Lepton Decay 

While the total two-photon decay rate provides in- 
formation about the constant part of the form fac- 
tors, the four-lepton rate can be utilized to probe both 
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FIG. 1: Double Dalitz Tree-Level Diagram. 
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FIG. 2: Double Dalitz Exchange Diagram. 



the momentum dependence of the form factors and the 
mixing of the couplings. The decay to four leptons 
may contain either two pairs of identical particles (e.g., 
Kl —* e+e~e+e~) or pairs of non-identical particles 
(e.g., Kl — » e+e^/i+/i^). The tree-level Feynman di- 
agram for the decay is depicted in Fig. |^. If the final 
state does contain identical particles, then there is also 
an exchange diagram, as shown in Fig. ^j. The matrix ele- 
ment is then the sum of the two diagrams A4 = A^i + A^2 
and its square is |A^p = |7Wi p -h |7W2p + 2 Re (A^| A^2)- 
The analysis presented in the rest of this section applies 
only to the direct contribution to the double Dalitz pro- 
cess. Appendix |C| gives an explicit expression for the 
interference term. 

The matrix element for the direct contribution to the 
double Dalitz process can be written as 



(12) 



where i? is a two-photon coupling given by Eq. (BS), k 



is a photon momentum, 11 is the propagator for a photon 
of momentum k 



(13) 



q and spin r and positron of momentum p and spin s 

r'^(g,p; r, s) = -ieUr{q)Yvs{p)- (14) 

The sum in the propagator extends over the three helicity 
states and the term proportional to fc'^fc" vanishes when 
contracted with the current. The matrix element can 
then be cast as 



A21 



(15) 



Ai Aq 



where Ti. is one of the two-photon couplings given in 
Eqs. (B4,B5). The quantity L contains all the lepton 



information and is equal to Lx = ie'^T ^/k"^ . 

The squared matrix element, summed over final state 
helicities can be expressed in terms of the five phase 
space variables a:;i2, X34, yi2, 2/34, and (j) (defined in Ap- 
pendix Focusing on the dependence, the squared 
amplitude is 



287r2a^ 



{A SVC? 



C sin (j) cos 4> + D sin ( 
£;cos0 4- F), 



(16) 



where 



A ^w^fp cos2 CX^[1 + (1 - A?2 + yl,)il - Xj, + y^J] 
+ fs sin' C^'[2 - - + Vu + 2/34]}, (17a) 

B =w^fl sin' C^'[l + (1 - A?2 + y?2)(l - AI4 + 2/34)] 
+ fp cos2 CA'[2 - Xj, - Xl, + yl^ + yj,]}, (17b) 

C =2fpfs sin ( cos ( cos SXzw'^ 

X (A?2 -y?2)(A34 -2/14) 

D =2fpfs sin ( cos ( cos SXw^yiy2 



yi2)(A34 - 2/34)' 



(17c) 



(17d) 



and r is the fermion current for an electron of momentum 



E =2/1 sin' C^^«' 2/122/34 V'(A22 - 2;f2)(Ai4 - 2/14), (17e) 

^^=/|sin^C^4(l-y?2)(l-2/i4), (17f) 

where the kinematic variables w, z, and A are defined 
in Appendix ^ Terms A and B arise from the diago- 
nal pieces of the transverse part of the pseudoscalar and 
scalar couplings while term F is due to the diagonal piece 
of the longitudinal part of the scalar coupling. Term C 
is the interference between the transverse parts of the 
pseudoscalar and scalar couplings while terms D and E 
are due to interference between the longitudinal part of 
the scalar and the transverse parts of the pseudoscalar 
and scalar couplings respectively. 

The partial width is obtained by integrating 1/(2M) 
times the square of the matrix element over the eight 
dimensional phase space (given in Appendix ^). The in- 
tegrals over yij extend from —Ay to +Xij and can be 
done analytically. The differential partial width, normal- 
ized to the two-photon width and integrated over j/12 and 
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TABLE I: Results of numerical integrations over X12 and 3:34 



assuming 


',fp = .fs = l 








Integral 


Kl — > eej-Ljj, 


Kl fJ-tJ-fJ-tJ- 


Kl — ^ eeee 


TT — > eeee 






9 Q77 . 1 




7 9987 


h 


0.11147 


1.149 ■ lO"-"^ 


14.201 


7.2838 


h 


0.21713 


1.120 • 10"^ 


28.343 


14.509 


h 


0.74203 


5.499 • 10"'' 


27.725 


15.600 


h 


0.76948 


1.595 • 10"^ 


27.809 


15.684 


h 


0.01503 


5.246 • 10"'' 


0.0556 


0.0555 



TABLE IL Coefficients of (j> dependencies for various values 
of the mixing angle (5 = 0). 





Kl — > ee/i/.i 


Kl fJ-fJ-fJ-fJ- 


Kl — > eeee 


tt" eeee 


RiO) 


0.7952 


1.379 • 10"'' 


8.6995 


4.8037 




0.8402 


5.313 • 10"'' 


8.7412 


4.8453 


Kl(0) 


-0.0668 


-0.0027 


-0.2033 


-0.1881 


Ki(V2) 


+0.0663 


+0.0027 


+0.2031 


+0.1879 


Co 


44.32° 


26.98° 


44.94° 


44.89° 




+0.0664 


+0.0026 


+0.2031 


+0.1880 



2/34, IS 



1 dVi Sa^ 



T^j d(j) Stt^ 



[h cos^ C sin^ (t> + h sin^ C cos^ 
+ I3 sin C cos C cos S sin cos 4> 



(18) 



where 5 is a symmetry factor which is 1/4 for modes with 
identical particles and 1 otherwise. The li factors rep- 
resent the integrals over X12 and X34 given below. The 
factors Ii and I4 correspond to the pseudoscalar cou- 
pling, I2 and I5 are the analogous terms for the scalar 
coupling, Iq is the additional longitudinal term in the 
scalar coupling, and I3 is the interference term. 



h = 



2 '^12'^34^^ 



•^12 '^34'^^^ 



^3 =g / / dxi2dx34fpfs- ^2 



(19a) 
(19b) 
(19c) 



Ii = J J dxi2dx3ifj> 

h = JJ dxi2dx3ifs 
1 



2 ^12A34A'^ , 2 



w 

2 ^12^34 



2 (3-A^2-A2j, (19d) 



-{3-XI2-XI,), (19e) 



h^^JJ dxi2da;34/iAi2A34A(3 - Af2)(3 - A^^). (19f) 

The double integral is performed by first integrating over 
X34 from 0:34 to (1 — ^2:12)^ and then over X12 from 



'12 

to (1— •\/x34)^, where = [mi + mjY /M"^ . In order to 
obtain values for these integrals, fp and fs must first be 
specified and then the integrals can be done numerically. 
Table | summarizes the values for the different double 
Dalitz modes assuming that fp — fs = ^- 

The numerical value of I5 + /g was found to be several 
orders of magnitude larger in Ref . . Extracting le from 
that result yields a value of 3578.0, compared with our 
value of 0.01503. This discrepancy has been traced to the 
use in Ref. ^ of a scalar coupling g^pgua, rather than 
the correctly antisymmetrized tensor g^pg^a — dfiagi^p- 
One consequence of the much smaller value we find for 
Jg is that the total width for the double Dalitz decay 



is almost completely insensitive to the mix of scalar and 
pseudoscalar couplings (except for the Kl /i+/i"/i+/i" 
decay), in contradistinction to the conclusion of Ref. j|] 
but in agreement with the comments near the end of 
Ref. §. 

The differential rate can also be expressed in a compact 
form, suitable for experimental fits to the 4> distribution, 
involving a constant term, a CP conserving cos 2(f) term, 
and a CP violating sin 2(j} term 



1 dTi _ 



i?(l + KiCOs2(?!) + K2sin20), (20) 



where 



R =5[(/i/2 + /4)cos2C 

+ (/2/2 + /5+/6)sin2C], (21a) 

Kl =S{l2 sin^ C - ^1 cos^ C,)/{2R), (21b) 

K2 =5/3 sin C cos C cos (5/(2 i?). (21c) 

The values of R and ki at C = and C — 7r/2, along with 
the maximum value of K2 and the angle ^0 at which it 
takes on that value, are listed in Table ||. As expected, 
for a pure pseudoscalar decay, the amplitude of the cos 2<^ 
term will be negative while the amplitude of the sin 2</) 
vanishes. For a pure scalar decay, the amplitude of the 
cos 2(j) term is nearly the same magnitude as in the pseu- 
doscalar decay but positive, and the amplitude of sin 2(j) 
again vanishes. Depending on the mode, the amplitude 
for the CP-violating sin 2(f) term is maximal for values of 
the mixing angle between 7r/8 and 7r/4. 

Alternatively, we could have integrated over before 
X12 and 2:34, in which case we would have 



1 d^Ti 2Sa^ A12A34A 



(3-A?2)(3-A24) 



^-y^ dxi2dx34 9tt'^ w'^ 

X [fp cos^ CA^ + /I sin^ C(A2 + 3wV2)], (22) 



where we have used — \^ 



The interference term 



integrates to zero and what remains clearly shows the 
kinematic differences between the contributions of the 
two couplings. 

Going back to Eq. (20), the final integral over (f) from 
to 2tt can be performed to get the direct contribution 
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TABLE III: The decay rate for pseudoscalar couplings — 0) 
assuming fp — fs = 1. 



Mode 


ri+2 




ri2/r^^ 




Kl — > ee/x/i 


2.859 


■ 10"^ 







2.859 ■ 10" 


-6 


Kl MA'A'A' 


9.914 


^0-10 


-0.512 • 10" 


10 


9.402 ■ 10" 


10 


Kl — * eeee 


6.256 


• 10"^ 


-0.036 • 10" 


-5 


6.220 ■ 10" 


-5 


tt" —^ eeee 


3.456 


• 10"^ 


-0.036 ■ 10" 


- 5 


3.420 • 10" 


-5 



to the four-lepton decay rate relative to the two-photon 
rate 



r 



2a^ 
37r2 



■R. 



(23) 



The total tree-level rate can now be computed for ar- 
bitrary form factors for modes without identical particles 
in the final state. For the other modes, there is the in- 
terference between the direct and exchange graphs that 
must be included. The decay rate has the form 



r = ri + r2 



(24) 



where, for modes without identical particles T2 — Vyi = 
0, and for modes with identical particles T2 = Ti. The 
expression given in Appendix ^ for the interference term 
could in principle be integrated numerically. We choose 
instead to use a Monte Carlo (MC) simulation to inte- 
grate the rate and make histograms of the relevant phase 
space variables. The decay rates for the various modes, 
broken into diagonal and interference terms, are listed in 
Table |nl 

Ref. 1^ included a similar table of values, some of which 
are in disagreement with our results. The most signifi- 
cant discrepancy involves the size of the interference term 
for the decays Kl — » e+e^e^e^ and 7r° —f e+e^e+e^. 
We find that the interference in Kl e^e~e~^e~ is 
roughly 9 times smaller than Ref. |^ reports, and that the 
interference in tt*^ — > e'^e~e^e~ is about 4 times smaller. 
We also differ in the total rate for Kl e"'"e~/z"'"/z~, but 
the factor of 2 difference is likely due to a typographical 
error in the previous publication. 

The assumption that the form factor is fiat contra- 
dicts current experimental findings. The two models that 
have been used to parameterize the kaon form factor are 
the BMS model || and the DIP model The BMS 
model was originally proposed to describe the coupling in 
the single Dalitz decay Kl e"'"e~7 and was therefore 
written as a function of one only 



1 



1 - rpX 



1 - rK'X 
1 

~9 



1 



1 - TpX 



1 



1 



1 



(25) 



The quantities r,; = M'^/Mf for Af,; equal to the p, K* , 
tu, or (j) meson masses. To apply this model to the double 



Dalitz decay, it is assumed that the coupling factors so 
that f{xi,X2) = f{xi,0) ■ f{x2,0)- In this paper, we 
will use a simplified form of the DIP form factor, which 
involves only the p meson and two parameters 



f{xi,X2) = 1 + anip 
+ I3dip 



Xl 



X2 



Xl - M2 /M2 X2 - M2 /AP 

XIX2 



(xi - M2/M2)(x2 - M^/A'P) ■ 



(26) 



As will be seen in Appendix the BMS model can be 
expressed as a generalized DIP model involving the p, uj, 
and 4> vector mesons. 

Experimentally, the form factor has traditionally been 
linearized in the case of the pion with just a slope pa- 
rameter measured, while for the kaon, the BMS model 
has been used and values of a^* quoted. The conversion 
to the DIP parameters is easily done, using the world 
average for the kaon we will use anip — —1.5 and 
for the pion, aoip = —1.0. There is as yet no experi- 
mental sensitivity to Pdip and so we will use Pdip — 0. 
The effect of using the DIP model with these values of 
a is that the n'^ e+e^e+e^ rate increases by less 
than 0.4%, the Kl e+e^e+e^ rate increases by 6.5%, 
the Kl — * e"'"e^/i^/i^ rate increases by 56%, and the 
Kl — > /i"'"/U^/i+/i^ rate increases by 68%. It is clear that 
the assumption of a flat form factor is completely invalid 
for modes containing muons. 



III. HIGHER ORDER PROCESSES 

The tree-level double Dalitz process is 0{a^) since 
it contains two electromagnetic vertices. Higher order 
contributions to the double Dalitz rate contain one or 
more internal loops. There are three types of graphs 
that contribute at O(e^): the vacuum polarization, the 
vertex correction, and the 5-point diagram. A represen- 
tative diagram from each of these processes is displayed 
in Figs, and |[ respectively. There are two graphs 
for both the vacuum polarization and the vertex correc- 
tion, one for each pair, plus four graphs for the 5-point 
function. If there are identical particles in the final state, 
there are exchange diagrams and the number of graphs 
doubles. The interference between the tree-level diagram 
and the one-loop diagrams is 0{a^) and therefore con- 
tributes to the first order radiative correction to the dou- 
ble Dalitz rate. 

Both the vertex correction and the 5-point graph con- 
tain IR divergences, that is, divergences in the limit that 
the exchanged photon energy goes to zero. In order to 
handle this singular behavior, one must also consider the 
radiative double Dalitz decay X — > Z^Zj~Z^Z^7, in which 
one of the leptons internally radiates a photon. There 
are two contributions to this process, shown in Figs. ^ 
and ^. The radiative process diverges in the opposite 
manner from the one-loop graphs making the combined 
decay rate finite. 
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The combined process will be indicated hy X 
I2 I2 il) , where the radiated photon may or may not 
be detectable. The distinction between the non-radiative 
and radiative decays is an experimental issue and ulti- 
mately related to the hardware. We will use the energy 
of the radiated photon in the CM frame to differenti- 
ate between events with a hard photon, > i?cut, and 



events without, < En 



The cutoff is chosen such 



that photons with energies below the cutoff can have no 
significant effect on the 4-lepton acceptance. The rate 
for radiative events with soft photons will be added to 
the rate for non-radiative decays. This contribution is 
also 0{a^) and therefore must be considered along with 
the one-loop corrections. 

The double Dalitz differential rate to second order can 
therefore be expressed as 



(27) 



where Jbrom is the brenisstrahlung contribution due to 
radiative decays with photons below the photon energy 
cutoff and Svin is the virtual correction due to the inter- 
ference between the tree-level and one-loop diagrams. 
The virtual correction can be further decomposed into 
the contributions from the three one-loop diagrams 



= 5vp + Sv 



(28) 



where d^p is the correction from the vacuum polarization 
diagrams, 6vc is the correction from the vertex correc- 
tion diagrams, and S^p is the correction from the 5-point 
diagrams. 



IV. RADIATIVE DECAYS 

The radiative double Dalitz decay will only be consid- 
ered at tree-level. It is straightforward but tedious to 
write down the expression for the rate. The two contri- 
butions to the rate are shown in Figs. ^ and ^. For each 
process, there exists three additional diagrams where the 
photon is radiated off of the other leptons, plus four ex- 
change graphs if applicable. 

Our results for the radiative decay rates use a MC sim- 
ulation in which we calculate the amplitudes for each he- 
licity state using explicit representations of the spinors 
and polarization vector. An photon energy cutoff of 
400 keV in the CM frame is used for kaon decays while for 
pions, a cutoff of 100 keV is used. It is useful to define the 
quantity = m'l^/M'^, where m^e is the reconstructed 
four-lepton invariant mass, to distinguish between the 
radiative and non-radiative processes. In terms of this 
variable, the cutoff for both kaon and pion decays is at 
0.9985. Fig. H shows the distribution of X4e for 



^cut 

— > e^e~e~^e~ and Kl e^e~e^e~^ events. The 
large peak at X4e = 1 is due to non-radiative events. The 
part of the distribution which falls away from the peak at 
1 is due to radiated photons from the process of Fig. ^. 
The rising part of the distribution near x^e = is due to 




0.; 



0.4 0.6 



0.8 1 



FIG. 3: for Kl e e e e and Kl e e e e 7 
events. 



TABLE IV; Tree-level rates for radiative decays including 
both all radiation and only hard radiation such that x^e < 
0.95, with fp = fs = 1. 



Mode 




r4i7(a;4e < 0.95)/r77 


tt" — > eeeej 
Kl eeee'y 
Kl eefifij 
Kl fJ-fJ-fJ-fJ-'J 


6.614(1) ■ 10-** 
1.540(1) ■ 10^5 
3.279(3) ■ lO"'^ 
5.634(3) ■ 10-12 


2.055(1) ■ lO"*^ 
0.504(1) ■ 10-5 
0.873(1) ■ 10-^ 
0.346(1) ■ 10-12 



hard Dalitz photons from the process of Fig. |[ Lost due 
to bin size is the low energy cutoff at E^ = 400 keV and 
the high energy cutoff at x^e = 16ml/M^. 



Tabic IV shows the tree-level radiative decay rates for 
the four modes, with no form factor. The rates in the first 
column include photons of all energies, while the rates in 
the second column include only photons with energies 
large enough that x^^e < 0.95. This value of x^e is chosen 
to closely match the resolution on the four-lepton mass 
in current experiments. 



V. RADIATIVE CORRECTIONS 

The next four sections will describe the different contri- 
butions to the radiative corrections to the double Dalitz 
differential rate. The first three contributions arc rela- 
tively straightforward to determine and we will therefore 
only summarize the relevant formulas. The last contribu- 
tion, coming from the 5-point diagram, is considerably 
more difficult to calculate. In particular, numerical in- 
stabilities plague the evaluation of the tensor 5-point in- 
tegrals involving light leptons. The fourth section, along 
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with much of the appendix, will outline our procedure for 
obtaining this (usually) small but non-negligible contri- 
bution. We will present the full 5-point diagram correc- 
tions to the differential rate in closed form. In Ref. 
van Ncerven and Vermaseren reported a numerical in- 
tegral of the radiative corrections to the related two- 



photon process e+e 



'n" but did not present the 



corrections to the differential cross section, 
discuss in Sec 



VI 



As we will 

the effects of radiative corrections 
are much more important compared to form factor ef- 
fects when considering double Dalitz decays, because the 
values in the accessible phase space are much smaller 
than those typically probed in two-photon resonance for- 
mation with final state lepton tags. 



X 




FIG. 4: Radiative Diagram 1. 



A. Bremsstrahlung Correction 

The contribution to the double Dalitz differential decay 
rate due to the soft bremsstrahlung part of the radiative 
decay is defined as 



Shi-cmixi2, X34,yi2,y34:, <j>) = 



(29) 



where d^rbrcm/rf^'I' is the differential decay rate for the 
soft part of the radiative decay integrated over the pho- 
ton momentum with the constraint Ej < Ecnt- The full 
differential rate is 



1 



^^rb.e. = ^El>^-J^rf'^<i>(^^, (30) 



where is the four-body phase space differential and 
A^brcm is the matrix element for the soft bremsstrahlung 
contribution. If the photon energy cutoff is taken small 
enough, the matrix element can be approximated as 



My 



P2 ■ e 



Pi 



Pi ■ e 



P3 



P2 ■ k Pi- k pi ■ k 



P3 



Mt 



(31) 

where e and k are the radiated photon's polarization and 
momentum 4-vectors, respectively. There is one contri- 
bution from each of the radiative diagrams represented 
by Fig. H The other type of radiative process (Fig. ||) 
does not contribute in this limit. If the cutoff is small 
enough, the lepton momenta can be held fixed while the 
photon momentum is integrated out, with the result 



Jo 



cPk 



(32) 



X 




FIG. 5: Radiative Diagram 2. 



where 



P2 



Pi 



Pi 



Ps 



P2 ■ k Pi- k pi ■ k P3 ■ k 



, (33) 





2pi • P2 




2P3 • Pi 






2pi ■ Pi 




{Pl 


■ k)(p2 ■ 




{P3 • fc)(P4 


k)^ 


(pi 


■ k)(pi 


k) 




2p2 ■ P3 




2p2 ■ Pi 






2pi • P3 




{P2 


■ k)(p3 ■ 


k) 


{P2 ■k){pi 


k) 


(pi 


■ k)(p3 


k) 




Pl 




Pl 


Pl 




Pl 





{pi ■ fc)2 {p2 ■ fc)2 (p3 • fc)2 {pi ■ fc)2 ■ 

The correction can be expressed in terms of a sum of 
ten integrals which can be done in closed form 

^brem = 47ra [2/(pi , P2 ) + 2I{p3,Pi) + 2I{pi,Pi) 
+ 2Iip2,P3) - 2I{p2,Pi) - 2/(pi,p3) 
- I{Pl,Pl) - I{P2,P2) - I{P3,P3) - I{Pi,P4)], (34) 



where 



I{ki, kj 



Scut 



d^k 



{2T:f2Ek ih-k){kj-k)- 



(35) 



Each integral yields both a finite part and an IR diver- 
gent part which goes as ln(2i?cut/A) where A is the pho- 
ton mass which will be taken to zero after the divergent 
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terms are canceled against each other. The first two di- 
vergent terms wih be seen to cancel the divergent parts 
of the vertex correction, the next four cancel divergences 
in the 5-point functions, and the last four cancel the elec- 
tron self-energy divergences (which are included in the 
renormalized vertex function). For ki ^ kj 



I{ki,kj] 



In 



In 



f2E, 



4 



V A 



Lio 1 - 



Lio 1 



T 0+ 



Li2 1 - 



Lio 1 



T fl^ 

T 

Xij Xij 



where 



nf = 

T — 



(1 + ^i,jkl ± Kjkl)/'^, 

(1 + SjAkl ± Xj.ikl)/ C^CTij), 

O-y (1 + Sijkl) - (1 + <>],ikl), 

(Zy+Ay)/(1-Zy +(5y), 



(36) 



(37a) 
(37b) 
(37c) 
(37d) 



and the various (5, z, and A symbols are defined in Ap- 
pendix For ki = kj 

1 

4^ 



I{ki, ki) 



In 



2£'ci; 



1 

2A" 



In 



1 + 
1 - A,, 



, (38) 



where A^ is again defined in Appendix 

It will be enlightening to extract the IR divergent part 
of the soft brem contribution and express it in a way that 
will make the cancellation obvious. Collecting terms, one 
has 



2a 



1 - 



Zl2 



1 



2Ai2 



In 



12 



h^^^ln 

TT Al3 

-—In 

TT Aon 



2A34 
/ Zl3 + Al3 
\Zl3 - Ai3 
2^23 + A23 



In 



Z12 + A 
Z12 — A12 

2^34 + A34 

Z34 - A, 



34 



Z23 



A 



23 



a zi4 
- -—In 

TT Al4 

a Z24 , 
-—In 

TT A24 



^^14 + Ai4 
Zi4 — Ai4 

^24 + A24 
^24 — A24 



(39) 



As will be seen shortly, the first two terms cancel the 
divergent part of the vertex correction while the last four 
cancel the divergent part of the 5-point correction. 



It 



X 



It 



FIG. 6; Vacuum Polarization Diagram. 

then the squared matrix element to 0{a^) is 

2Re {Ml,,M^My 



1 



\Mt 



= |A^trccP(l +<5virt). 

This defines Svwt, 

2Re(M*,,,Xvirt) 



^virt, — 



\Mt 



(41) 
(42) 

(43) 



Therefore, we must compute the matrix element for each 
of the one-loop contributions. 



1. Vacuum Polarization 

The vacuum polarization process involves higher order 
corrections to the photon propagator and is a function 
of the square of the photon momentum, or the x of that 
pair. There is one contribution for each photon propa- 
gator. One contribution is shown in Fig. |^. The vac- 
uum polarization diagram is IR finite but UV divergent. 
The divergence can be handled by renormalization of the 
photon wavefunction. The vacuum polarization matrix 
element can be written as the tree-level matrix element 
times the renormalized polarization insertion 



A^vp = A^troc X! n/(a;y ), 



(44) 



where the sum is over lepton species in the loop and the 
renormalized polarization insertion is 



B. Virtual Correction 

As mentioned above, the virtual correction arises from 
the interference between the tree-level and one-loop di- 
agrams. If the full matrix element is 

M=M tree + M virt 

+ 0(e6), (40) 



^i{xij) = — / dz z{l- z) 
Jo 

X In [1 - z(l - z)xijM'^/mf - ie] , (45) 

where mi is the mass of the lepton in the loop. The inte- 
gration depends on the size of Xij compared to mf /M"^, 
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X 




X 




FIG. 7; Vertex Correction Diagram. 



FIG. 8: 5-Point Diagram. 



such that 



for Xij > imf/M'^, while 



In 



1 + 



a f 



37r [3 

Pij 



{3 + pI) [7r-2tan-V^j] 



(46) 



(47) 



for Xij < imf /AP. The functions p and /? are related to 
A as defined in Appendix ^ but are functions of the loop 
mass 



= .yi-4m2/(a:,,Af2), 
= ^Amf/ix,jM^) - 1. 
The correction then is 



^Reni^(xy) 



(48a) 
(48b) 

(49) 



where the first sum is over the number of vacuum polar- 
ization graphs and the second sum is over the possible 
lepton species in the loop. 



divergences will be handled simultaneously by renormal- 
ization of the electromagnetic coupling and the lepton 
wavefunction while the IR divergent part will cancel the 
IR divergence in the soft bremsstrahlung correction. The 
matrix element corresponding to one of the diagrams is 



^ij J Uij ) 



where 



V(a;,y,yy) = Fi{xij) + ^2(0;^) 



Fi and F2 are vertex form factors defined by 

„ , N ".11, — , I -11 I A) 

Fi{Xij) 



(50) 



(51) 




where Zij and Xij are defined in Appendix ^ and A is the 
photon mass. 

The correction is then just 



2^ReV(a;,j,y„) 



(54) 



2. Vertex Function 

The vertex function involves higher order corrections 
to the QED vertex and is a function of the momenta of 
the pair. One contribution is shown in Fig. The 
vertex correction contains both UV and IR divergences. 
We will also include the self-energy correction to the lep- 
ton lines which also are UV and IR divergent. Both UV 



where the sum over g is over the number of vertex func- 
tion graphs. 

The IR divergent part is contained within Fi. Twice 
the real part of the In A term is exactly what is necessary 
to cancel the first two divergent terms (one for each ver- 
tex correction graph) in the soft brem contribution. The 
final four divergent soft brem terms will have to wait for 
the last piece of the virtual correction, the 5-point dia- 
gram. 
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3. 5-Point Diagram an exchange process. The diagram shown in Fig. g con- 

tains a photon exchanged between and l'^ . The matrix 
There are four distinct 5-point diagrams that con- element for that graph is 
tribute to the direct process, plus four more if there is 



dH 2 

(2^ M ^■^p^^"'P'' + ^sigt^pg^a - gp<ygyp)] {pi2 + tTipzi - ty 



{pi2 + t)2 - A2 + ie J |_ (p34 - ty ~ X^+ie 



i2 - A2 -I- ie 

X u{p2){-iejs) i 77-7 — — I {-ie-fa)v{pi)u{p4){~iejp) ( - — - — — ) {~iej,^)v{p3), (55) 

where t is the foop momentum and pi, p2, Pa, and p4 are the momenta of l^, , I2 , and I2 , respectively. This can 
be re-expressed as 

^s{[{Pi2 ■ P^dg'^" ~ P^^P-^] [hoA^. + I^iB^.a + I^ic^^ap] + {p34 - Pi2)p.g^. [l^^Af^" + 1^^ B^, + I^fc^p] 
- gp. [{h2r,A^'' + ihsY^pB'^'"' + {h4tp,C^^'"'] }) , (56) 



r 



where A, B, and C are combinations of spinors and 
gamma matrices and p^ — Pi + P2 + Pa + Pi- The fac- 
tors of li are integrals over the loop momentum. There 
are three basic integral forms from which all the other 
may be obtained. The notation for the integrals has the 
following meaning: the first digit in the subscript refers 
to the number of denominators and the second refers the 
number of powers of the loop momentum appearing in 
the integral. The 5-point integrals are defined in Ap- 
pendix ^ as a function of four 4-vector arguments fci, 
k2, ks, k^. 

For the diagram shown in Fig. the arguments of the 
5-point integral functions should take on the values 



ki = -P2, 



k2 = ~Pi, 
ki = -pi, 



so we can write 



ho{-P2,-Pl,P5,-P4:) 



(57a) 
(57b) 



(58) 



for the scalar integral, and analogous expressions for the 
higher-rank tensor integrals. The spinor terms for this 
diagram are 

= -4(p2 • P3)u{p2h^v{pi)u{p4)Yvip3), (59a) 
B'"-" = 2[uip2)rv{pi)u{pi)Y-/''MP3) 

- uip2)hl'^^''vipi)u{Pi)Yv{p3)], (59b) 

C^'^P- = u(p2)7''7^7^«(Pi)«(P4)7"7"7^'*^(P3). (59c) 



The spinor terms for the diagram containing a photon 
exchanged between and I2 are 

= 4(p2 • P4)u{p2h''v{pi)u{p4h''v{p3), (60a) 
B^-P = 2[u{p2)^ij^Yv{p,)u{pi)Yv{p3) 

- u{p2)j''v{pi)u{p4)^27'^YviP3)], (60b) 

C^-P- = -u{p2)j''n'v{Pl)u{pihr,Yl^<P3), (60c) 

and the scalar integral for this diagram is 

ho = l5(){~P2,~Pl,P5,~P3)- (61) 

The spinor terms for the diagram containing a photon 
exchanged between if and I2 are 

= 4(pi • P3Mp2h''v{pi)u{pih''v{p3), (62a) 
B^'^P = 2[u{p2)YY'MPiMP4h''viP3) 

- u{p2)^''v{pi)u{pi)Yl^i>iv{p3)\, (62b) 
Ci^upa ^ -u{p2)Yri''v{pi)u{pih^^''j^v{p3), (62c) 

and the scalar integral for this diagram is 

ho ^ l5o{-Pi,-P2,P5,-P4)- (63) 

The spinor terms for the diagram containing a photon 
exchanged between if and I2 are 

= ~4(pi • pi)u{p2h^'v{pi)u{p^)j''v{p3), (64a) 

B^-P = 2[uip2)rv{pi)u{Pi)^iY'l'<P3) 

- u{p2)Yl'^i>iv{pi)u{pi)Yv{p3)], (64b) 
CP'-P- = u{p2)Yri''v{pi)u{pi)^^-i''Yv{p3), (64c) 
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and the scalar integral for this diagram is 
ho ^ 4o(-Pi, -P2,P5, -Ps)- 



(65) 



The tensors A, B, and C are computed for a given 
helicity combination and combined with the integrals to 
yield the matrix element for that helicity state. The cor- 
rection then involves a sum over the sixteen possible final 
states 



Eiii2Re[A^r„e(A)A^5p(A)] 



35p 



where A here refers to the helicity state and 

A^trcc(A)=^XLc(A), 

a 

A^5p(A)-5]A^fp(A). 



(66) 

(67) 
(68) 



The sums here are over the number of graphs for each 
process. 

The IR divergent part of the 5-point correction is most 
easily isolated by looking at the 5-point matrix element 
in the IR limit. All terms involving tensor integrals van- 
ish leaving only the /50 term. The divergent part of /50 is 
due to the two divergent box integrals, J^q^ and I^^ . The 
relevant terms in the scalar 5-point function for Fig. ||, 
in terms of the divergent 3-point function, are 



rlR _ _ - 
^50 - 2 



1 



4j 



Mi 



P?2 



Af2 



= W[(P?2-Mf)(P34-M|)] • (69) 
Extracting the In A piece of ///j, II^q can be written as 



tIR _ 
-'50 — 



16Tr^X23PhiPi2 
In 



'Mf)ipl, 
Z23 + A23 



^23 



A2 



Ml) 



- 2iTX 



in A, 



(70) 



where /i is a kinematic function with dimensions of mass 
which is independent of A. The divergent part of the 5- 
point matrix element, dropping the finite term involving 
is proportional to the tree-level matrix element 



a Z23 



•27r A 



23 





+ A23\ 




— A23 / 




\Z23 





In A. 



(71) 

The IR divergent part of the 5-point correction coming 
from all four diagrams is 



'Sij — ^iTL 



IT Xi 



A,; 



A,; 



In A, 



(72) 



where Sy is the product of the sign of the charges of pi 
and pj and the sum is over the four diagrams. Again, it 
can be seen that this is the necessary form to cancel the 
remaining four divergent brem terms. 




FIG. 9: Distribution of the total radiative correction for tt 
eeee events with an IR cutoff of 3:4"* = 0.9985, with fp 



TABLE V: Average size of radiative correction to the differen- 
tial rate with an IR cutoff of xlf = 0.9985, with /p = /s = 1- 





tt" eeee 


Kl — > eeee 


K]^ eefifi 




s 


-0.1948 


-0.2618 


-0.0788 


+0.0805 



VI. MC SIMULATION RESULTS 

The inclusion of the radiative corrections impacts both 
the differential rate and the total rate. The total correc- 
tion to the differential tt^ ~> e^e~e^e~ rate is shown in 
Fig. ^. The average size of the correction factor for the 
four different modes is shown in Table The total rate 
for the combined 4-lepton plus photon process is inde- 
pendent of IR cutoff. Table VI summarizes the tree-level 

cutoff independent 



rate and the rate for the combined, 
process, divided into the rate including all radiation and 
the rate including only soft radiation (x4e > 0.95), all 
with fp = fs — ^- It is the last column which should 
most accurately predict the observed non-radiative 4- 
lepton rate. It is seen that the non-radiative rate is 
smaller than the tree-level rate for both 4-electron modes 
while it is larger for the modes with muons. 

The probability of radiation can now be computed as 
the ratio of the radiative rate to the combined rate. 



Table VII lists the probability of radiating a photon 
(x4e < 0.9985) along with the probability of radiating a 
hard photon {x4e < 0.95) for each of the four modes. The 
probability is highest for — > e^e~e''"e~ (7) where the 
X values can be the smallest. The probabilities for Kl 
e'^ e~ p,~ {■y) are slightly less than half of what they are 
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TABLE VI: Summary of tree-level 4-lepton rate and com- 
bined radiatively corrected 4-lepton plus photon rate, includ- 
ing the rate for all X4e and the rate for x^^ > 0.95, using 
fp = fs = 1. 



Mode 


T^tree /-p 
^41 1 ^ 11 


^ii{i)/^ii 


all X4,e 


a;4e > 0.95 


tt" eeee 
Kl eeee 
Kl eefifi 
Kl fJ-l^l^l^ 


3.421(4) ■ 10"^ 
6.222(5) ■ 10~^ 
2.858(1) • 10'*' 
0.941(1) ■ 10"^ 


3.536(4) ■ 10"^ 
6.406(4) ■ IQ-^ 
2.996(3) ■ 10"^ 
1.026(1) ■ 10"^ 


3.331(4) ■ 10"^ 
5.903(4) ■ 10"^ 
2.909(3) • 10"^ 
1.025(1) ■ lO"'' 



TABLE VII: Probability of radiation {xi^ < 0.9985) and 
probability of hard radiation [x^e < 0.95), defined as P = 

Yil-1 1^4,1(1) ■ 



Mode 


P{x4^ < 0.9985) 


P{x4e < 0.95) 


tt" eeee(7) 


0.187 


0.058 


Kl eeee{'j) 


0.240 


0.079 


Kl eefifi{j) 


0.109 


0.029 


Kl ^ MMMM(7) 


0.006 


0.0003 



M (7), 



for the four electron mode, and in 
there is very httle radiation. 

The effect of the radiative corrections on the differ- 
ential rate can be observed in the distributions of the 
five phase space variables. The statistics in the follow- 
ing plots reflect the amount of CPU time dedicated to 
each mode. While the calculation of the radiative cor- 
rections is CPU-intensive, it is actually the generation 
of the radiative decays that takes the most time. 

For the modes with identical leptons, it is useful 
to adopt a method of pairing the electrons with the 
positrons in order to study the dilepton mass distribu- 
tions. We choose to use the pairing for which the product 
of x's is minimized. It is this pairing that will contribute 
the most to the matrix element in general. Therefore, 
Xa and Xh are the z's belonging to this pairing, with the 
additional requirement that Xa < Xb- In addition, ya is 
the y variable defined in the a-pair CM, and yb is the 
same quantity in the 6-pair CM. And lastly, (f>ab is the 
angle between the planes of the a-pair and &-pair in the 
overall CM. 

The first variable that we will look at is x which is 
modified by both the existence of a form factor and 
the inclusion of the radiative corrections. In all cases 
we set /3 = in the DIP form factor model. Figs, 
and |ll| show the distribution of Xa and Xb, respectively, 
for 7r° ^ e+e~e+e~ events. The plot on the left com- 
pares the distribution using the tree-level matrix element 
with no form factor (a = 0) to that using the same ma- 
trix element but with a — —1.0. The plot on the right 
compares the distribution using the tree-level matrix el- 
ement with no form factor to that using the radiatively- 
corrected matrix element also with no form factor. We 
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FIG. 10: (a) Xa in vr ~-^ eeee events using the tree-level differ- 
ential rate with q = (dots) and with a = —1.0 (histogram), 
(b) Xa in tt'' eeee events using the tree-level differential 
rate with a = (dots) and the corrected rate for events with 
X4e > 0.95, also with a = (histogram). The ratio of the 
dots to the histogram in both cases are shown in (c) and (d). 



have provided a linear fit to the ratio over some reason- 
able range on a scale appropriate for comparing the two 
effects. For the form factor comparisons, the dependence 
should be primarily linear. This is not the case for the 
radiative corrections in general. The P^r degree of 
freedom is included as a measure of the linearity. It can 
be seen that the form factor has a much smaller effect 
on the X distribution than the radiative corrections do. 
This is not too surprising since the range of accessible 
values for the 7r° decay is relatively small in addition to 
being far from our assumed p pole. 

For the kaon modes, we observe that the form factor 
has a much larger effect on the x distribution than the 
radiative corrections do. Figs. |l2| and |l^ show the distri- 
bution of Xa and Xb, respectively, for — > e+e^e+e^ 
events. The plot on the left compares the distribution 
using the tree-level matrix element with no form factor 
(a = 0) to that using the same matrix element but with 
a —1.5. The plot on the right compares the distri- 
bution using the tree-level matrix element with no form 
factor to that using the radiatively-corrected matrix el- 
ement also with no form factor. The roll off at high x 
in plot (c) of Fig. |l^ is due to presence of the exchange 
diagram in this mode 



Figs 

tributions for Kl — > e^e^/i ' /i 
no pairing ambiguities and we plot Xee and a;^^. It can 
be seen that there is no roll off in plot (c) of Fig. |l^, and 
furthermore, a small quadratic dependence is observable. 
While the x of the ee pair is slightly modified by the ra- 
diative corrections, the x of the fj,fi pair does not change 
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and y_5| show the same dis- 
events. Here there are 
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FIG. 11: (a) Xb in tt" — > eeee events using the tree-level differ- 
ential rate with a = (dots) and with a = —1.0 (histogram). 

(b) Xb in tt" eeee events using the troe-levcl differential 
rate with a — (dots) and the corrected rate for events with 
X4e > 0.95, also with a = (histogram). The ratio of the 
dots to the histogram in both cases are shown in (c) and (d). 



FIG. 13: (a) xi, in Kl — > eeee events using the tree-level 
differential rate with a = (dots) and with a = —1.5 (his- 
togram), (b) Xb in Kl eeee events using the tree-level 
differential rate with a = (dots) and the corrected rate for 
events with X4e > 0.95, also with a — (histogram). The 
ratio of the dots to the histogram in both cases are shown in 
(c) and (d). 
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FIG, 12: (a) Xa in Kl eeee events using the tree-level 
differential rate with a = (dots) and with a = —1.5 (his- 
togram), (b) Xa in Kl — > eeee events using the tree-level 
differential rate with a = (dots) and the corrected rate for 
events with X4e > 0.95, also with a = (histogram). The 
ratio of the dots to the histogram in both cases are shown in 
(c) and (d). 



FIG. 14: (a) Xee in Kl ee/i/i events using the tree level 
differential rate with a = (dots) and with a = —1.5 (his- 
togram), (b) Xee in Kl — » ee/i/i events using the tree-level 
differential rate with a = (dots) and the corrected rate for 
events with X4e > 0.95, also with a = (histogram). The 
ratio of the dots to the histogram in both cases are shown in 
(c) and (d). 
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FIG. 15: (a) x^fj, in Kl — > ee/x/x events using the tree-level 
differential rate with a — (dots) and with a — —1.5 (his- 
togram), (b) Xfj_fi in Kl eefifi events using the tree-level 
differential rate with a = (dots) and the corrected rate for 
events with X4,e > 0.95, also with q = (histogram). The 
ratio of the dots to the histogram in both cases are shown in 
(c) and (d). 



FIG. 16: (a) ya in Kl — * eeee events using the tree-level 
differential rate (dots) and the corrected rate for events with 
X4e > 0.95 (histogram), (b) yt for the same events. The ratio 
of the corrected distributions to the tree-level distributions 
are shown in (c) and (d). 



at all. This is as expected for the massive muons. 

Fig. |l^ shows the distribution of ya and yi, for the tree 
level differential rate and the radiatively corrected dif- 
ferential rate for Kl e+e^e+e^ events. The effect 
here is quite small. Since y is a measure of the energy 
asymmetry of the lepton pair, it is seen that the radiative 
corrections tend to make the pairs slightly more asym- 
metric on average. 

The effect on the cj) distribution is due entirely to the 
5-point diagram. Fig. ^ shows a comparison of the dis- 
tribution of (pab generated with the tree-level matrix el- 
ement to the same distribution generated with the ra- 
diative corrections, for Kl I-J,'^ ^j.^ ij.'^ fi^ events. The 
enhancement at (j)ab — tt and the corresponding depletion 
at (pab = = 27r can be understood in terms of Coulomb 
interaction between the final state particles. The config- 
uration at (/) = has all leptons in a plane with the same 
sign particles near each other. The effect is only observ- 
able in the Kl — *■ ^l^ ^~ decay where the leptons 
in each pair are usually well separated. 
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FIG. 17: (a) The distribution of 4>a,b in Kl MMMM events 
using the tree— level differential rate (dots) and the corrected 
rate for events with x^^ > 0.95 (histogram), (b) The ratio of 
the corrected distribution to the tree-level distribution. 



VII. CONCLUSIONS 

The main conclusion that can be drawn from these 
distributions is that the radiative corrections are ex- 
tremely important for extracting a form factor in the 
tt" — !■ e"*"e^e+e~ mode. For the kaon modes, the form 
factor has a larger impact on the x distribution and 



the modification of the distribution due to the radia- 
tive corrections is less important. The only published 
result for the Kl e"*'e~e"'"e~ mode Q quotes aoip = 
— 1.1 ± 0.6(stat). Fig. ^ shows that at most the radia- 
tive corrections would change the slope by 0.22, which 
while significant, is smaller than the current experimen- 
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tal error. Likewise, for the Kl e+e~/x+/x~ mode, 
the latest result ||l^ based on the invariant mass shape is 
ctDiP — 53^2 70- T^^^ present experimental error is 
again larger than the impact of the radiative corrections 
on the mass distribution. 

As for the extraction of the mixing angle C from the 
observed (j) distribution, the radiative corrections can be 
safely neglected at present. 

The two publications above quote an integrated rate 
(normalized to the two-photon rate) of (6.24±0.34) • 10~^ 
for Kl -> e+e-e+e- and (4.51 ± 0.42) • 10"^ for Kl -> 
e^e~ (jf^ li~ , where the errors are purely statistical. These 
results are in good agreement with our predictions when 
both the radiative corrections and a form factor with 



-1.5 are included. For Kl 



e ' e e ' e 



, the two 



effects offset and the net result is an increase of just less 
than 2% over the tree-level rate with no form factor. In 
Kl — > e~ fi^ fj,~ , the form factor is the dominant effect. 
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APPENDIX A: KINEMATICS 

The four particle final state can be kinematically de- 
scribed by considering subsystems containing only two 
particles. Consider the system composed of two par- 
ticles with momenta pi and pj and total momentum 

pfy We will 



Pi + Pj and mass squared 



Ptj 

define a dimensionless dot product of any two vectors pi 
and Pj as 



where 



Zij ^2{p,-pj)/plj = l-x 



Xt^pf/pij-, Xj^p^pi^. 



(Al) 



(A2) 



The energy and momentum of each particle in the two- 
particle CM frame are 

E* = my (1 + hj)l% E* = rn,, (1 - 6,,)/2, (A3) 
p*=myAy/2, (A4) 

where 



Wi 



2^a 



(A5) 
(A6) 
(A7) 



Occasionally we will use symbols like Zijk whose meaning 
is interpreted as z{pi,pj +pk)- 



Now consider a three body system composed of mo- 
menta Pi, Pj, and pk- There are two phase space vari- 
ables needed to describe the system. The first one will 
De Xij — Pij /pjj^. The other one is defined in the ij CM 
frame as the cosine of the angle between the direction of 
particle i and particle k. 



cos 9. 



ij,k 



"^Pk ■ (Pt -Pj) 



' Pijk^ij^ij-k 



Pijk^ij^ij,k 



(A8) 



A more convenient variable that will be used in place of 



COS! 



:,fc IS 



= Xij COS I 



ij.ki 



(A9) 



where k will always refer to the total momentum minus 
the ij momentum. 

Finally, the four body final state requires five phase 
space variables to uniquely describe it. We will use the 
X and y values for the two lepton pairs plus the angle 
between the normals of the planes defined by each pair 
in the overall CM frame. The first four variables are 



=p?2/A/^ 

X3i^pljM^, 

2p34 • (pi - P2) 



yi2 



2/34 



APX 

2pi2 • {P3 - P4) 



(AlO) 
(All) 

(A12) 
(A13) 



where the second term in the numerator of the j/'s van- 
ishes since S12 = <^34 = 0. Any use of z, w, or A without 
subscripts will refer to the functions of X12 and X34, so 
A = Ai2,34 for instance. The last phase space variable is 
defined as 



(j) = ta,n ^ {smcj}/ coscp). 



where 



sm ( 



cost 



I6€f,i,pap'^pi^pf^p1 



M^Xw^{Xi2-y-U)i>^U-yLy 
M'^zyuysA - 2(pi - P2) ■ (ps - P4) 
M^w^{Xl2-yl2){Xl,^yl^) 



(A14) 

(A15) 
(A16) 



The angle </> is defined so that at = the two pairs lie 
in a plane with the like-sign particles adjacent to each 
other. The orientation of ^ = tt again has both pairs in 
a plane, but with the opposite signed particles adjacent. 
The general expression for the phase space integral is 



1 (Ppi(Pp2d^P3d^P4 



fO \8 TtiT? 5^{P~Pl~P2-P3~Pi)- 

(27r)» l&Ep^Ep^Ep^Ep^ 

(A17) 

Upon integrating out the (5-functions, integrating over 
the Euler angles, and changing variables to those listed 
above, the phase space reduces to 



214^ 



Xdxi2dx34dyi2dy: 



'34 1 



(A18) 



16 



where the factor 5 is a symmetry factor which is required 
for modes containing identical particles in the final state. 
The double Dalitz modes with identical particles contain 
two sets, thus requiring two factors of 1/2. So S — 1/4 
if the final state contains identical particles, and S — 1 
otherwise. 

When there are identical leptons in the final state, the 
amplitudes for the exchange diagrams have the same al- 
gebraic form as for the non-exchange diagrams except 
that the kinematic variables Xi2, 0:34, ?/i2, yai, and (j) 
are replaced by 2:14, X23, Uu, 2/23. and ^14,23- These ex- 
change variables will in general be functions of all five 
of the non-exchange variables. As is turns out, we only 
need explicit representations for 2:14 and 0:23. These are 
given by 

xiA = [1 - Ai22:i2 - A34.X34 - A(yi2 - yu) - zyuyu 



and one finds only two couplings 



+ (A?2 - yh){^li - yli) cos (/.]/4, (A19) 
a;23 = [1 - Ai2a;i2 - A34X34 + \{yi2 - ysi) - ^2/122/34 

+ [Xh - 2/?2)(Ai4 - yh) cos^]/A. (A20) 



APPENDIX B: MESON-77 COUPLINGS 

In this section we will work out the explicit form of the 
two-photon couplings, allowing for photons of arbitrary 
mass, using the polarization vectors in the helicity basis. 
The general form of the coupling is 



'H\i\2 — ^A'i'pCT^i'eAi^2^A2 ' 

where H is either 



M 



P ^fiupaj 



(Bl) 

(B2a) 
(B2b) 



The three polarization vectors for a massive photon in 
the helicity basis are chosen to be 

e+(±i) = (0, l,±i,0)/V2, (B3a) 
e-(±i) = (0,l,Ti,0)/V2, (B3b) 
e°(±£) = (fc, 0, 0, ±E) /%/fc2, (B3c) 

for a photon traveling in the ±£ direction. With these po- 
larization vectors, one finds three couplings for the scalar 



-MTsz ,Ai=A2 = +, 
^fiA. = <! -MTsz , Ai = A2 = -, (B4) 
-MTsw ,Ai=A2 = 0, 

where z and w are defined in Appendix |A[ The longi- 
tudinal contribution vanishes for the pseudoscalar case. 



-iMTpX 
-iMTpX 







, Ai — A2 — +, 
, Ai = A2 = — , 
,Ai = A2 = 0, 



(B5) 



where A is also defined in Appendix |^. There are three 
interesting differences between the scalar and the pseu- 
doscalar couplings. First, assuming that (5 = 0, there is 
a relative phase between them. Additionally, the trans- 
verse couplings differ in the kinematic factor. And lastly, 
there is the additional scalar coupling due to the contri- 
bution from longitudinally polarized photons. Where as 
the transverse couplings go like A or z, both of which are 
0(1) on average, the longitudinal coupling goes hke w, 
which is 0{x), making its contribution less significant. 



APPENDIX C: DOUBLE DALITZ 
INTERFERENCE 

The interference between the tree-level direct and ex- 
change contributions for modes with identical leptons is 
a sum of three terms 



2Re(7W*A^2) 



M'^Xi2X'iAXlAX2:i 

X {Afp{xi2,x:ii)fp{xi4:,X2z)cos^C, 



+ B[fp{xi2,x:ii)fs{xii,X2-i) 

+ /p(a;i4,a;23)/5(a;i2, 2:^34)] sinCcosC 

+ C'/s(a;i2,X34)/s(a;i4,a;23)sin2Q, (CI) 



where 



A = A2{2S2 + 877^ + S [8772 - z(yi2 + y34)'] 
~w'^{\ + yi2y34){2 - 2/12 - 2/34) 
+ ■^'>f{yi2 + y3A){xi2yi2 + 2^342/34)}, (C2) 



B = -2A[ry2(a;i2 -|- a;34 - Xu - ^23) 

-I- 'i{xi2X3i - a;i4a;23)]Stan(/), 



(C3) 



C = -~~?z - ~?[&Tfz -w^- (2z2 + W^)yi2y34] 
~ E[z{z^ - 2w^) + z{z^ + 2w^)yl^yli + 2z^yi2y34 
- 2Tf{z'^ + w^){l + 3J/122/34) + 877^2] 

I 2 2 3 3 22 I 2/0 2 r, 2\ 2 2 

+ w z 2/122/34 -wz yi22/34 + w (3z - 2w )2/i22/34 

- 2zi;2(z2 _ ^2)(y2^ ^ y2^^ ^ ^2(^2 _ ^^2^ 
r, 2ro 2 2 2 0/ 2 3 2\ 

- 277 [3w 2:?7i22/34 - 2(z ~ z -w )yi22/34 

+ 2z{z^ - 3wV2) - 2{z^ + w^)] 
-8ij%z-z^)-w^yi2y34], (C4) 



17 



where 2 = ^^/(Afa - 2/^2) (^^34 " vidcoscj) and 77^ = 
Aw? /M'^. The exchange variables a; 14 and 0:23 are de- 
fined in Appendix ^ in terms of the five non-exchange 
phase space variables. The term proportional to cos^ C, re- 
sults from interference between a pseudoscalar coupling 
in both the direct and exchange graphs, while the one 
proportional to sin^ C is due to scalar couplings in both 
graphs, and the one proportional to sin C cos ^ is due to a 
pseudoscalar coupling in one graph and a scalar coupling 
in the other. 



nominator. One can express tensor integrals in terms 
of lower rank tensor integrals with the same number of 
propagators and lower rank tensors with fewer propaga- 
tors |23|| . In the end, every tensor integral can be decom- 
posed into scalar 2-, 3-, 4-, and 5-point functions. The 
scalar 5-point function is not independent and can itself 
be expressed in terms of scalar 4-point functions. 



APPENDIX D: 5-POINT FUNCTION 

The matrix element for the 5-point diagram is com- 
posed of tensor integrals with five propagators in the de- 



This appendix will outline our procedure for first re- 
ducing the tensor integrals to scalar integrals, and then 
computing the scalar integrals in closed form. 



a. Tensor Integrals 

Begin by defining 

f dH 1 f dH f^t" 

[ d'^t [ d^t t^'ftP 

I^,{ki,k2,k,,k,)^J I,P,^P{k,,k2,k,,k,)^J (D2) 

where 

Ni^t^~fil N2 = {t + kxf - 111, N-i ^ {t + ki + k2f - til, (D3) 
Ni^ [t + ki+k2 + ks,f ~ ^xl, N5 ^ {t + ki+k2 + k3 + k4f - (D4) 

where fii is an internal mass and the ki are external momenta. 

The original reduction scheme of Ref. ]2^ , while theoretically sound, suffers from uncontrollable numerical inaccu- 
racies. To avoid this problem, we follow the procedure suggested in Ref. j2^, and use a reduction scheme based on 
the Schouten identity which utilizes Gram determinants to express any tensor integral as a sum of integrals, one with 
the same number of propagators and the rest with one less propagator, and all with the power of the loop momentum 
reduced by one. The identity has the following form 

^/.gfelfe2fe3fc4 ^ . fc^)gA'fefc3fe4 . J^^^^^ktt^k^ki _^ (-^ . ^3)gfelfe2Mfe4 . ) gfel fe2 fcs M ^ (D5) 

^/J _ 1 ^j<q_^^tJ.k2kiki _ jy^^tiiki+k2)k3k4, _|_ J^^^tiki{k2 + k3)k4, _ J^^^p,kik2{k3+k4,) ^ ^^g/jfci fe2 fes ^ (D6) 

where 

Z;^ = (S . fc^)eMfe2fe3fc4 _ . J^^yi^k^k^k^ ^ . ^^y^J.k^k2k^ _ . fc_j)£Mfelfe2fc3^ (QJ) 

and is defined in terms of its dot products 

■s ■ ^1 = ^ [m2 - Ml - ^i] : s ■ ^3 = ^ [m4 - Ml - (^1 + ^2 + kzY] - s ■ ki - s ■ k2, (D8) 

■s ■ ^2 = ^ [m3 ^ Ml - {ki + fe)^] - s - ki, s ■ k4 = ^[f4- fJ-i- {ki +^2 + ^3 + ^4)^] - s ■ ki - s ■ k2 - s ■ k^. (D9) 

The notation e'=i''2'=3fe4 jg shorthand for e''"'"' ki^,k2uk'ipkia- 
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The reduction then proceeds as follows: 
1 



^51 



^52 



^53 



ckik2k3ki 



1 



ckik2k3ki 



1 



cki k2 ks 



50 



Ufe2fe3fe4 7-(l) ,/j(fei+fc2)fe3fe4 r(2) i ('=2 + fe3 )fc4 /"(S) .M*: 1 '=2 (fcs +*:4 ) 7(4) 

t J40 t J40 -|- e e 



^fikik2k3 jiji) 



-t^kikaki I j-Wf i,vtW] Mki+k2)k3ki t{2)i' , nki(k2+k3)ki AS)^ 
' M4I ~'*^l-'40 41 -'41 

_ gMfeife2(fe3+*;4)jW'' _|_ ^tiktk2k3 j{5)v I 

nk2k3ki( j{'i-)vp l,v tWp _ l,P tW'^ , uuip A^)] Mkl+k2)k3ki t{'2)vp 
fc 1 J42 '4'1-'41 '^1-'41 "r'^l'^l-'40 I ^ -'42 

^ ^tiki(k2+k3)ki j^(i)'^P _ gMfeife2(fc3 + fe4) j(4)i'P _|_ ^nkik2k3 j{5)vp 



(DIO) 



(Dll) 



(D12) 



where for any ho = J (2^*4 n,N2N3N, ^ 



TP _ 
-'41 — 



9j;fcife2fe3 
^"kik2k3 



"fclfc2fc3"fclfe2fc3-'40 

M2 — 2jfclfe2fe3 1 felfc2fc3"felfe2fc3-'41 
/ci ^2 ^3 



<rp.k2k3 r(l) _ rMCsi +^2 ) ^3 7-(2) , rAlfel (fc2 +^3 ) 7-(3) _ r^lfel fe2 7-(4) 
"fclfc2fc3-'30 "feifcsfes 30 +"fcifc2fe3 -'30 "fcifc2/C3 30 



cnk2k3 I TV l.vA^)\ ^p.{ki+k2)k3 A2)v 

"felfc2fc3l 31 '^^l-'aO ) "fcl/C2/C3 -'31 



rM'=l('=2+fc3) 7-(3)l' _ cMfelfe2 7-(4)l' 
' "felfc2fc3 -'31 ''fcife2fe3 31 



2gfclfc2fe3Mgfelfe2fc3f 

efelfe2fc3 
felfe2fc3 



0sfcifc2fc3-'40 + 2 ' 



(s + fci)fe2fe3 t(1) _ rs(fci+fe2)A;3 7-(2) 
fclfc2fc3 -'30 "fcifc2fc3 30 



.sfcl(/C2+fc3) r(3) _ csfeifc2 r(4) 
"felfc2fe3 30 "feife2fc3 30 



(D13) 



(D14) 



and for any /30 = / 



t 1 

(2;0^ WiAr2Af3 



— 



fcifc2 



'felfc2"fcifc2-'30 ^ 



^^"^2 ,-(1) _ cA'(fel+fe2) r(2) cpfci ^(3) 

"fcife2-'20 "fcifcs 20 +"fcifc2-'20 



(D15) 



The reduction notation has the following meaning: /4Q is the 4-point function obtained from /50 by dropping the 

ith propagator, /3Q is the 3~point function obtained from its corresponding 4-point function by dropping the jth 
propagator, and so on. The Gram determinants that appear in the reduction are kinematic functions which are 
defined as 



skik2 _ 
9192 



qi ■ ki qi- k2 
92 ■ ki q2- k2 



<:kik2k3 _ 
9192 93 



qi ■ ki qi ■ k2 qi ■ k^ 
q2 ■ ki 52 ■ k2 92 • h 
qs ■ h 93 • k2 93 • ^3 



s:kik2k3k4 _ 
9192 93 94 



91 ■ ki 91 ■ k2 91 • ^3 91 • k4 

92 ■ ki 92 • k2 92 • fcs 92 • ^4 

93 ■ ki 93 • k2 93 • ^3 93 • ki 
qi ■ ki 94 • k2 94 • ^3 94 • ^4 



(D16) 



The traces that appear in Eq. pq can also be reduced 



{h2rp 
(43)7 



j(i) 
-'40 ' 

_ J,aAl) 
J4I t-i ^40 J 



r(l)«/3 
'42 



k" T 

"-1 



41 



^1 ^1 ^ ^ 



40 



(D17a) 
(D17b) 
(D17c) 



While this procedu re is generally much more reliable then that of Ref . , there are still problems that occur when 
as defined in Eq. (D7) is not an independent combination of the final state momenta. This happens when all the 
momenta in the parent particle's rest frame lie in a plane. Even though these configurations form a subspace of zero 
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P5 



\ i - Pi - P2 




FIG. 18: Scalar 5-Point Function. 



measure in the final state phase space, finite numerical precision dictates that they will be generated with non-zero 
probability by the MC program. In this case, we use the identity 



1 r 



where = g'=i'=2fc3M^ xhe first term on the right vanishes upon integration, allowing l!^-^ to be written as 



7-p 

^51 



9 rfcife2fe3 



rsa/3 cfj.a0 J 



clik2k3 t(1) _ eAtfe3fe4 r(2) ,(-^^4^1 7-(3) _ i;/ifcife2 r(4) 

"felfe2fe3 40 "feifc2fe3 40 "fclfe2fe3 40 "feife2fe3 40 



(D18) 



(D19) 



The other tensor integrals can be expanded in a similar manner. The same problem can arise in the reduction of the 
4-point tensor integrals if the three momenta ki, k2, and ks are linearly dependent, in which case this same procedure 
is reproduced at one lower rank. 

In these degenerate cases we have a choice between numerical inaccuracies resulting from antisymmetric invariants 
such as e'=i'=2fe3fe4 being very small or inaccuracies resulting from assuming exact linear dependence. To decide which 
approximation is be tter, we do the calculation of the tensor integrals both ways and check whether identities such as 
those in Eqs. ( D17 ) are satisfied. More than 99% of the time, one of the two methods yields good agreement for all 
these "trace checks" . 



b. Scalar Integrals 



The most general 5~point function we will need to consider is 



ha = Ai 



1 



d^t 

(27r)^ 1 [f2 - fj,j + ie][{t + fci)2 + ie][{t + ki + kaf - f4 + le] 

1 

^ [(t + ki+k2 + ksf - m1 + *e] [{t + ki+k2 + k3 + k^f - ^lI+k 



(D20) 



In the case at hand /ii = 0, = Wi, /is — Mi, fi4 — M2, and /is = m2 where the m's are lepton masses and the Af's 
are boson masses. In addition, ki = —p2, ^2 = —pi, ^3 = Ps, and k^ = —p4, where the external momenta satisfy the 
relation pi + P2 + Ps + Pi = P5- The diagram representing this integral is shown in Fig. 18, 



We have allowed the boson propagators to have arbitrary masses Mi and M2 in order to include a form factor in the 
calculation of the 5-point function. This is necessary since the form factor becomes a function of the loop momentum. 
We will use a generalized DIP form factor which, with the appropriate choice of coefficients, can reproduce both the 
DIP and the BMS form factor models. The generalized DIP form factor is 



/(a;i,a;2) 



1 



Xl 



X2 



xi-M^/M^ X2-Mf/M^ 



X1X2 



{xi - M^/M^){x2 



M^/M^ 
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where M is the total mass and the sum is over propagator masses Mi. The diagram containing two photon propagators 
and a form factor is then replaced by a sum of "n? diagrams of four different types, one containing two photon 
propagators, two containing one photon and one massive boson propagator, and one containing two massive boson 
propagators. 



X1X2 



1 



X1X2 



1 



1 



xi{x2-Mi/M^) X2{xi- M'i/M^) 



+ 



E 



A: 



{xi-Mf/M^)[x2~M^lM^Y 



(D22) 



The general form of Eq. ( D20 ) permits the evaluation of all four of these contributions. To use a specific model, the 
parameters and /3y must be adjusted. For the simple DIP model, only a p meson term is included with ai = a dip 
and Pii — [3dip- The DIP model requires the inclusion of four 5-point diagrams involving all combinations of photons 
and p mesons. The BMS model is more complicated, requiring 25 different diagrams. To simplify this we have let 
= Mp, which reduces the number of diagrams to 16. The values of the generalized parameters, in terms of aK* 

are 



Q!l = —1 



^{l-Mj/Ml,) 



0.2 



1 



9 \\-MllMl, 



Ml/Ml.' 



(D23a) 

(D23b) 

(D23c) 



and Pij — aiUj. A flat form factor is obtained by se tting all of the ai — Q and (3i.j — 0. 

We will write the 5-point function in Eq. (D20) as a sum of 4-point functions using the following relationship 
between n-point functions and (n — l)-point functions |23], 



Xc./«, + (n-5 + 2.)co/,f=«-^^ 

i=l 



(D24) 



where 



Co 



and 



0, 



•5.. =(a^?+M--4)/2. 



for i < j. 



(D25) 
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We will use to refer to the propagator mass when th e dist inction between vector bosons and leptons is irrelevant. 
In the case at hand n = 5 so the second term in Eq. (D24) is 0{e), and since the 5-point function in D = 6 — 2e 
dimensions is finite in the limit e — > 0, the scalar 5-point function can be written as a sum of five scalar 4-point 
functions 



1 ^ 

Je 



5 



7 = 1 



-1 



'40 



(D27) 



a. The 4^-Point Function Two of the five 4-point functions contain IR divergences due to the presence of the 

3- point functions where both of the vector boson propagators have been removed. Therefore there are two distinct 

4- point functions that we will need. The first has one zero mass propagator 



1 



{2tt)D [t2 + je][(t + kif - + + fci + ^2)^ - + i€][{t + ki + k2 + k^f -^ll + ie 



(D28) 



and the other has two zero mass propagators and two lepton propagators 
I'iQ{ki,k2,kz,mi,m,2) = 

,2c 



1 



{2Tr)D [i2 + + ki)^ + + k, + -ml + te][{t + fci + fcz + ^3)^ - + ie 
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where in I'^q, (fci + fc2)^ = mf and (fci + k2 + k^)"^ = m\. 

In order to use Eq. ( D28| ) for where all four propagators have non-zero masses, and to extract the divergent 



part of l'^ and l'^ , we make use of the following propagator identity 



1 — a 



[{q + pY - nW [{q + p + kf - f4] [{q + p)^ - ^il]{q + p + ak^ [{q + p + k)^ - fij] {q + p + ak)^ ' 
where a is chosen to be the positive root of the equation 

a{l - a)k^ - (1 - Oi)^\ - anl = 0. 
This identity allows us to write the five 4-point functions as 

= aiho{-aiP23, -pi,P5,mi,Mi,M2) + (1 - q;i)/4o((1 - ai)p23,P4 - P5,P5,m2, Mi, M2), 

(2) 

= ho{-Pi2,P5,-P4,Mi,M2,m2), 
I4Q = a34o(~"3P34, -P2,P23,'rni,m2) + (I - a3)^4o((l - a3)P34,Pl -P5,P23,M2,mi,m2), 

= a44o("4Pi2,P3, -P23,rn2,mi) + (I - q;4)/4o(-(1 - a4)pi2,Pi,P23, A'A, Jtii, TO2), 
I40 = ho{-P2, ~'Pi,P5,mi, Ml, M2), 



where 



ai = (l + 523 + A23)/2, 



as = I - MI/pI^, 



a4 = 1 - M2/p?2 



(D30) 

(D31) 

(D32a) 
(D32b) 
(D32c) 
{D32d) 
(D32e) 
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The finite 4-point function /40 defined in Eq. ( D28 ) can be expressed in closed form as a sum of 36 dilogarithms. 
We will define it in terms of the function 



J (A, B) = / [ln(z -B±ie)- \n{A -B± ie)] . 

Jq z-A 



= Lio 



A 



A- B±ie 



Li? 



A-1 



A- B±ie 



(D34) 



For arbitrary complex arguments, A and B, the integration would also produce additional logarithms with prefactors 
which depend on the relative difference between the signs of the imaginary parts of A and B ||2^ , however if A is real 
these additional terms vanish. In the case at hand, A will always be real and we will only need the dilogarithms. In 
terms of these new functions. 



ho{kl,k2,k3,^il,^i2,^J'3) 



1677^7(77+ — 77 ) 

— e — d + ie^ 



J 



rj- 



-d + ie 



1-/3' e + k 



J 



J[(3 + V 
-r/+ —d + ie 



-e-d 



k 
J 



J 



yy"*" —d + ie 



1-/3' e + k 



J {f3 + r]+,z+ +ie) - J {l3 + j]+,z^ - ie) + J {(3 + r] , z+ + ie) + J {l3 + rj , z^ - ie) 



J 
J 



T^'"2 



zj + ie] +J 



+ -J 



1-/3' 2 



Z2 - if ] - J 
J 



V 



1-/^2 



1-/3' 



,z^ + ie\+J 



If 



1^3 



13 ' ^7 " V /3 ' ^7 ' " V /3 

where P is either root of the equation + jP + b ~ 0, rj^ are the roots of the equation 

[eg — jk — I3gk]rj^ + [eh — ck — dj — 2(3dg]rj + [ae — cd + j3{ak — dh)] — 

and ^ — eg — jk — Pgk. The quantities zf^ are the roots of the following equations 

= gzl + {h + j + k)zi + {a + b + c + d + e), 
0= {b + g + j)zl + {c + e + h + k)z2 + [a + d), 
= bzi + {c + e)z3 + (a + d). 



(D35) 

(D36) 

(D37a) 
(D37b) 
(D37c) 
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The lower case variables are combinations of the elements of the relevant 4x4 matrix Sij defined in Eq. ( D25| ), 



a — 5*33 + 5*44 — 25*34, 

b = S22 + S33 — 2523, 

C = 2(523 — 524 — 533 + 534), 

d = 2(534 — 544), 
e = 2(524 - 534), 



/ = Si4, 

g — 5ii + 522 — 25i2, 

h = 2(5i3 — 5i4 — 523 + 524), 

j = 2(5i2 - 5i3 - 522 + 523), 

k = 2(5i4 - 524). 



(D38a) 
(D38b) 
(D38c) 
(D38d) 
(D38e) 



The divergent 4-point function I'^q defined in Eq. ( D29 ) can also be written in closed from. The divergent part is 
just the divergent 3-point function 



lIR 



dH 



(27r)4 [t2 - A2 + ie] [{t -p2f-m\ + ie] [{t + pj.f-ml + ie] ' 

^23 + A23 A ^. ^ , vlt . o . ( 2A2 



327r2p23A 



23 



In 



4«7r In A23 — In 



, 2i7r 

^23 — A23 

1 + <523 ~ A23 



ln^~27r2 + 2Li2 



+ 2Li2 



2A 



In^ 

1 
2 



1 ^ ^23 + A23 



In 



2A23 

2 / 1 + (^23 + A23 



23 



In^ 



In 



1 + <^23 + ^2 



1 + ^23 + A23 

1 — <523 ~ A23 



2A 



23 



In 



1 + <523 + A23 



2A- 



23 



In^ 



1 — <^23 + A23 



In A23 In 



-In^ 



2A23 

1 ~ ^23 + A23 

2A^^ 

1 + (523 ~ A23 



1 — ^2Z + A23 
1 — '^23 + A23 



In 



In A23 In 



23 ^23 



-In^ 



2A23 

1 — <523 ^ A23 
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The full expression is then 



-f4o(^l>^2,fc3,"^l,"^2) = 



+ Li2 



167r2A23A:| 
\-A 
B- A + ie 



In 



2fci • ^3 



fc2 



In 



- Li2 



1-A 



C-A-ie 



Z23 — A23 

^23 + A23 

- Li2 



-f 2z7r 

-A 



B - A + 



+ Li2 



-A 



C-A-ie 



ln(l - A - ie) [ln(l - B- ie) - ln(l - C + ie) - \n{A - B - ie) + \n{A -C + ie)] 
\n{-A - ie) [\n{-B - ie) - ln(-C + ie) - ln(^ - B - ie) + ln{A -C + ie)] \ ] , (D40) 



where 



A^-{kl/2 + h-k2)/{h-k3), 



S = (f + (523 + A23)/2, 



C = (f + ^23 - A23)/2. 



(D41) 



It is worthwhile to consider the special case of two photon propagators, that is Mi = M2 — 0. In that case the 
4-point functions simplify considerably. We will write 



4o(^l>^2,/C3,AiliM2) 



dH 



1 



(2^)4 [i2 _ ^2 ^ ie] [{t + fci)2 - /i2 ^ [(^ ^ ^ ^2)2 + ie] [(i + fci + ^2 + A:3)2 + ie] ' 



Li2 

hLi2 

Li2 

HLi2 



V 



+ 



- Li2 



1 



- Lio 



— — ze 



V — ri 

- 



— Li9 



v+ -1 
w+ — ri 

Li2 

' — V 

V- -1 



w+ — r2 



- Lio 



V — — le 



V — 7'1 

Li2 



-Li2 



- Li2 

V- - 1 

- 7'2 
+ Li2 



- Lio 



^v^ — v 
,+ - 1 
— + ie 

V 



V — V 

- 1 



le 



where A is the discriminant and are the roots of the quadratic equation 

[k{h + k) - dg]v'^ + [e{h + k) - dj]v -bd^O, 



(D42) 
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and are the roots of 



gv'^ + {h + j + k)v + b = 0, 



and 



ri — {—d—e + ie)/k, 
Therefore, when Mi = M2 — 0, the five 4-point functions are simply 



^2 = ie/{h + k). 



(D44) 



(D45) 



j(l) 

-'40 


= 4o(P23,P4,-P5,"ll,"i2), 


(D46a) 


,-(2) 
-'40 


= 4o(P3,P4,-P5,0,TO2), 


(D46b) 


7-(3) 
-'40 


= -^4o(~-P34,-P2,P23,'71i,m2), 


(D46c) 


j(4) 
-^40 


= Iio{Pi2,P3,-P23,fn2,mi), 


(D46d) 


j(5) 
-'40 


= 4o(--P2,-Pl,P5,0,TOi). 


(D46e) 



Also, in this case, /4Q sim plifie s somewhat because A defined in Eq. (D41) becomes one. The first two dilogarithms 
in the second line of Eq. ( D40 ), along with the entire third line, vanish in this case. I'^q then becomes 



lio{ki,k2,k3,mi,m2) = ^( Iir + ttt^. — 7^ 

Kf \ 167r^A23«:3 



— Li9 



-1 



In 

+ Li2 



2fci • fcg 



ln(^Hi^|+2^7^ 



2:23 + A 



23 



-1 



B-l + ieJ \C-l-ie 
+ ill [ln(-B - ie) - ln(-C + ie) - ln(l - B - ie) + ln(l - C + ie)] 



(D47) 



where B and C are still given by Eq. (D41). 

b. The 3-Point Function There are ten 3-point functions that are needed, all of which are finite except one, 
/jq''^ = Iin defined in Eq. ( D39| ). The superscripts used in this section denote the two propagators that have been 
dropped from the original 5-point function to obtain the particular 3-point function. The finite 3-point functions 
can be generically written as 



-f30(fcl,fc2,Ml>M2,Ai3) 



1 



{2tt)D [i2 _ ^2 ^ + + te][{t + fci + fc2)2 - f4 + le] ' 



(D48) 



167r2(c + 26/3) 
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+- U2 
+ Li2 

+ Li2 

— Li2 



vi — v-^ + ie 
f 2 - 1 
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V3 

— v'^ — ie 
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where the Vj are 



Vl 



2a + d + l3{c + e) 
c+26/3 ' 



V2 



(l-/3)(c + 26/3)' 



V3 



V3 — + ie 



d + e(3 
(3{c + 26/3) ' 



and the vf^ are roots of the three quadratic equations 



0^bil + (c + e)ii + {a + d + f), 
= (a + 6 + c)i>'i + {d + e)v2 + f, 
= avl + dii3 + /, 
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and /? is either root of the equation 6/3^ + c/3 + a = 0. The lower case letters are again combination of the relevant 
matrix Sa 



a — S22 + S33 — 2S23, 
b = Sii + S22 — 2S'i2, 



C = 2{Si2 — Si3 — S22 + S23), 
d = 2(523 - S33), 



e = 2(513 - 'S'23), 
/ = ^'33. 



The ten 3-point functions that we need can then be expressed as 



c. 



7-(l2) 
-'30 


= ho{P4,-P5,m2,M2,Mi), 


-^30^^ = hoiPs, -Ph + Pi, 0, m2, Ml), 


/■(13) 
-'30 


= ho{P23,P4,mi,m2,M2), 


4o'^ =-r3o(-Pi2,P5,0,Mi,M2), 


r(14) 
-'30 


= ho{P23,-P5 +P4.,mi,m2,Mi), 


7-(34) _T^^ 
^30 — ^IR, 


7(15) 
-'30 


= ho{-Pi,P5,mi,Mi,M2), 


= ho{-P2,P5 - Pi, 0, mi, M2), 


x(23) 
-'30 


= ho{P3,P4,0,m2,M2), 


= l3o{-P2, -Pi, 0, mi, Ml). 


The 2-Point 


Function And finally, the general expression for the 2-point function is 



l20{kl,Hl,H2) 
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Tie) 



2 
1 

47r^' 



(27r)^ [t2 - ^2 + + fc^-)2 _ ^2 + ' 
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+ 2 - (1 - u+) In (1 - v+ - ie) - (1 - v') In (1 - t;" + ie) 



— \n {—v~^ — ie) — V ln{—v + ie) 
where u^are roots to the quadratic equation 

ky + (^2 _ ^2 _ ^2)^ + ^l=0. 
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(D53a) 
(D53b) 
(D53c) 
(D53d) 
(D53e) 
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The UV divergent term containing e cancels when the 2-point functions are combined to form the tensor integrals 
and can therefore be safely ignored. The ten 2-point functions that we require are then 
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-'20 


= /20(-P4,M2,TO2), 
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-'20 


= /2o(-pi,mi,Mi), 
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= ^2o(P3,0,m2), 
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= /20(P34,0,M2), 
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= /20(-f^l2,0,Mi), 
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-'20 


= l2o{P5- Pl,mi,M2), 
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= -^2n(-P2- 0. mi). 


(D57c) 
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